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Cohesive  Zone  Approach  to  Multiscale  Modeling  of 
Nanotube  Reinforced  Composites 

Executive  Summary 

As  a  part  of  the  project,  a  number  of  scientific  advances  were  made.  Though  this  project 
covered  funding  only  one  year,  it  provided  the  necessary  background  material  for 
advancing  the  state  of  the  art  on  cohesive  zone  model.  One  of  the  difficulties  in  cohesive 
zone  model  is  that  the  model  parameters  cannot  be  accurately  evaluated  using 
experiments.  Molecular  dynamics  provided  that  data. 

ABSTRACT 

In  order  to  fully  harness  the  outstanding  mechanical  properties  of  carbon  nanotubes 
(CNT)  as  fiber  reinforcements,  it  is  essential  to  understand  the  nature  of  load  transfer  in 
the  fiber  matrix  interfacial  region  of  CNT  based  composites.  With  controlled 
experimentation  on  nanoscale  interfaces  far-off,  molecular  dynamics  (MD)  is  evolving  as 
the  primary  method  to  model  these  systems  and  processes.  While  MD  is  capable  of 
simulating  atomistic  behavior  in  a  deterministic  manner,  the  extremely  small  length  and 
time  scales  modeled  by  MD  necessitate  multiscale  approaches.  To  study  the  atomic  scale 
interface  effects  on  composite  behavior,  we  herein  develop  a  hierarchical  multiscale 
methodology  linking  molecular  dynamics  and  Unite  element  method  through  atomically 
informed  cohesive  zone  model  parameters  to  represent  interfaces.  Motivated  by  the 
successful  application  of  pullout  tests  in  conventional  composites,  we  simulate  fiber 
pullout  tests  of  carbon  nanotubes  in  a  given  matrix  using  MD.  The  results  of  the  pullout 
simulations  are  then  used  to  evaluate  cohesive  zone  model  parameters.  These  cohesive 
zone  models  (CZM)  are  then  used  in  a  finite  element  setting  to  study  the  macroscopic 
mechanical  response  of  the  composites.  Thus,  the  method  suggested  explicitly  accounts 
for  the  behavior  of  nanoscale  interfaces  existing  between  the  matrix  and  CNT.  The 


developed  methodology  is  used  to  study  the  effect  of  interface  strength  on  stiffness  of  the 
CNT  based  composite. 

1.  Introduction 

Carbon  nanotubes  possess  excellent  combination  of  mechanical  and  thermal 
properties  with  high  aspect  ratio.  Stiffness  of  CNTs  measured  experimentally  1 1 1  and 
calculated  from  simulations  [2]  is  of  the  order  of  1000  GPa,  while  the  nearest  competitive 
fiber  (SiC  whiskers)  has  utmost  400  GPa  stiffness  [3],  CNTs  have  tensile  strength  of 
about  150  GPa  [4],  The  range  of  elastic  deformation  and  fracture  strain  are  also 
extremely  high  [5,6].  In  addition,  the  typical  length  of  CNTs  is  of  the  order  of  a  few 
microns  while  the  diameters  range  from  less  than  nanometer  {for  single  wall  nanotubes) 
to  about  30  nm  (for  multi  wall  nanotubes)  [7,8]  that  correspond  to  an  aspect  ratio  around 
1000.  This  combination  of  mechanical  properties  has  raised  the  great  possibility  of 
obtaining  super-strong  and  stiff  composites  with  CNTs  as  reinforcements.  Further,  the 
excellent  electrical  and  optical  properties  of  CNTs  facilitate  development  of  many 
multifunctional  products  [9].  Bulk  composites  and  thin  films  of  CNT  based  composites 
with  a  variety  of  polymer  matrices  such  as  polystyrene,  PMMA,  Polyurethane,  PPA, 
epoxy  and  PmPv  are  being  fabricated  by  a  host  of  researchers  [8,  10-14],  In  addition, 
metallic  and  ceramic  materials  such  as  Cu,  Al,  Ti,  SiOj  and  Al;Oj  have  also  been  used  as 
matrices  [15-19]  with  CNTs  as  reinforcements. 

Good  properties  of  fibers  do  not  necessarily  translate  into  good  properties  of 
composites.  There  are  several  issues  pertaining  to  alignment,  dispersion  and  load  transfer 
which  have  to  be  optimized  in  order  to  achieve  the  possible  high  strength  and  stiffness  in 
the  composites.  This  is  especially  true  when  the  reinforcing  elements  are  in  the 


nanometer  dimension  as  in  CNTs.  it  is  well  established  through  a  number  of  years  of 
research  in  conventional  composites  that  interface  strength  is  a  critical  parameter  in 
controlling  the  load  transfer  between  fiber  and  matrix,  and  forms  the  key  issue  addressed 
in  this  work. 

There  are  conflicting  observations  in  the  existing  literature  regarding  the  nature  of 
interfacial  bonding  in  CNT  based  composites.  The  nanometer  length  scale  of  CNTs 
renders  controlled  experimentation  to  measure  the  interfacia!  strength  extremely  difficult, 
if  not  impossible.  As  a  result,  most  of  the  observations  have  been  extrapolations  based  on 
microslructural  characterization  and  bulk  property  measurement.  Several  researchers 
have  observed  evidence  for  load  transfer  based  on  microscopic  and  spectroscopic 
observations  [8,  10,  20,  21].  Wagner  and  coworkers  [11]  report  a  high  value  of  interface 
strength  of  500  MPa  based  on  fragmentation  tests  in  urethane  -  CNT  composites.  Cooper 
and  cow'orkcrs  [22]  report  variation  of  interface  strength  in  the  range  of  30-300  MPa 
based  on  CNT  pullout  using  scanning  probe  microscope.  On  the  contrary,  there  are 
reports  of  low  interfacial  load  transfer,  for  e.g.  Schaddler  and  coworkers  [23]  use  Raman 
spectroscopy  to  infer  low  load  transfer  between  concentric  lubes  in  multi  wall  nanotubes 
in  epoxy  matrix.  Similarly,  Ajayan  and  coworkers  [24]  observe  low  load  transfer  in  CNT 
ropes.  Most  importantly,  the  measured  values  of  mechanical  properties  of  CNT  based 
composites  reported  by  various  researchers  are  much  lower  than  those  predicted  by 
theoretical  models  [25].  For  example,  Andrews  and  coworkers  [12]  report  normalized 
Young’s  modulus  (stiffness  of  composite/stiffness  of  matrix)  values  in  the  range  of  1.2- 
2.5  for  polystyrene  matrix  reinforced  with  2.5-25  volume  %  CNT.  The  comparable 
theoretical  maxima  and  minima  based  on  parallel  and  series  models  are  in  the  range  of  14 


to  140  and  1.03  to  1.33  respectively.  Parallel  model  corresponds  to  a  perfectly  aligned 
reinforcement  with  perfect  bonding  and  thus  represents  true  upper  bound.  Series  model  is 
the  case  when  the  alignment  of  all  reinforcements  are  normal  to  the  load,  but  still  perfect 
bonding  and  represents  a  possible  lower  bound.  An  intermediate  value  is  obtained  if 
either  the  reinforcements  are  randomly  aligned  and/or  the  interfacial  bonding  is  less  than 
perfect.  As  the  experimental  values  are  closer  to  the  lower  bound,  one  of  the  likely 
reasons  is  the  lack  of  adequate  bond  strength  given  that  the  alignment  is  possibly  random. 
Thus,  we  have  a  practical  reason  to  understand  the  CNT-matrix  interfaces  in  these  classes 
of  composites. 

The  purpose  of  the  paper  is  two  fold.  First,  we  establish  the  fact  that  interface 
properties  at  nanoscale  can  be  significantly  improved  by  enhancing  the  bond  strength 
between  the  matrix  and  CNTs  by  chemical  attachments.  We  show  that  this  is  possible  by 
conducting  simulations  of  pullout  tests  on  a  single  CNT  fiber  embedded  in  a  matrix  using 
molecular  dynamics.  Second,  we  intend  to  show  that  these  enhanced  interfacial  properties 
at  the  atomic  level  translate  into  substantia]  increase  in  the  value  of  composite  stiffness  at 
the  macro  level.  This  entails  development  of  a  multi  scale  model  with  FEM  modeling  the 
macroscopic  continuum,  while  MD  simulates  the  nanoscale  interfaces  with  atomically 
informed  cohesive  zone  parameters  providing  the  necessary  coupling  between  MD  and 
FEM. 

It  should  be  noted  that  there  are  a  plethora  of  approaches  in  the  general  area  of 
multiscale  modeling  concerning  two  different  length  scales.  We  focus  our  attention  only 
on  cases  where  the  smaller  length  scale  is  at  atomistic  level  (described  by  atomic 
simulations)  and  the  larger  length  scale  is  at  macroscopic  level  (described  by  continuum 


concepts).  The  difficulty  encountered  in  this  class  of  problems  is  that  while  atomistic 
regime  is  conceived  as  discrete  set  of  atoms  (e,g.  modeled  by  molecular 
dynamics/statics),  continuum  methodologies  assume  the  presence  of  differentiable 
displacement  field  with  no  specific  assumptions  on  the  underlying  physical  matter  that 
causes  the  said  motion. 

Many  multiscale  methods  such  as  quasicontinuum  method  [26],  coarse  grained 
molecular  dynamics  (CGMD)  [27]  and  hand  shaking  methods  [28]  employ  intense 
atomic  calculations  near  the  core  of  atomic  scale  inhomogeneity  and  use  continuum- 
based  solutions  in  the  homogenously  deforming  regions.  The  two  regions  are  linked  by 
various  methods.  Quasicontinuum  invokes  Cauchy-Born  rule  in  the  continuum  regime 
while  CGMD  employs  statistical  averages  of  atomistic  variables.  In  the  handshaking 
methods  kinematic  restraints  are  applied  at  the  boundary  of  the  two  regimes. 

In  addition  to  above  multiscale  methods,  continuum  methods  have  been  applied  to 
study  CNTs  and  CNT  based  composites,  which  are  essentially  nanoscale  systems.  For 
example,  Yakobson  [5]  has  shown  that  predictions  of  continuum  shell  theories  are 
compatible  with  atomic  simulation  results  of  deformation  of  nanolubes.  There  have  been 
attempts  to  model  the  elastic  behavior  of  nanolubes  using  continuum  techniques  (e.g. 
[29]).  Li  and  Chou  [30]  model  nanotubes  as  atomic  trusses  with  stiffness  having  all  the 
components  of  covalent  bonds  (angle  dependence,  torsion  etc).  Wagner  and  coworkers 
[11]  have  used  Kelly-Tyson  model  and  experimental  results  to  obtain  interfacial  strength 
in  urethane  matrix  nanotube  composite,  Throntenson  and  coworkers  [31]  use 
micromechanistic  model  to  obtain  composite  effective  properties.  Fisher  and  co workers 


[32]  have  incorporated  nanotube  waviness  into  micromechanics  model  of  nanotube 
composites.  These  approaches  assume  perfect  bonding  at  CNT-Matrix  interface. 

The  difference  between  atomic  simulations  and  continuum  methods  is  with 
respect  to  both  length  and  time  scales.  While  the  elegance  and  advantages  of  unified 
multiscale  models  such  as  CGMD  and  quasicontinuum  are  undeniable,  their  application 
is  constrained  due  to  the  disparate  time  scales  of  the  underlying  processes  occurring  at  the 
different  length  scales.  One  approach  to  counter  this  problem  is  by  hierarchical  mulliscale 
modeling  in  which  the  information  of  atomic  simulations  is  processed  and  fed  into 
continuum  models.  In  this  paper,  such  a  hierarchical  multi  scale  model  to  address  the 
problem  of  interfaces  in  nanotube-based  composites  is  presented. 

We  use  molecular  dynamics  simulations  to  simulate  CNT  pullout  test  from 
polymer  matrix.  The  interface  is  modeled  as  hydrocarbon  chemical  attachments  bonded 
to  both  matrix  and  fiber.  These  simulations  are  used  to  obtain  the  constitutive  response  of 
the  interface  in  the  form  of  tract  ion -displacement  plot  for  cohesive  zone  model.  This 
atomically  informed  CZM  is  used  in  finite  element  framework  to  study  the  constitutive 
response  of  the  composite,  thus  traversing  across  the  length  scale  from  nano  to  macro 
scales. 

Because  of  the  hierarchical  approach,  we  can  examine  various  factors  at  atomic 
scale,  which  affect  the  form  and  magnitude  of  parameters  that  are  passed  on  to  continuum 
simulations.  Since  CZM  is  central  to  our  model,  section  2  explores  how  interfaces  are 
represented  in  both  continuum  and  atomic  scales.  Section  3  presents  the  atomic 
simulations  of  pullout  test  and  explains  the  method  of  extracting  the  CZM  properties 


from  MD  simulations.  Section  4  discusses  the  application  of  the  multiscale  model  to  the 
problem  of  determining  the  stiffness  of  the  composite. 

2.  Problem  Formulation 

Interface  in  general  can  be  defined  as  a  region  that  separates  two  distinct  bodies. 
The  bodies  in  question  may  be  of  same  material  (e.g.  crystals  with  different  orientation) 
or  different  materials  (as  in  the  case  of  composite  materials).  The  bodies  may  react 
differently  to  themiomechanical,  electromagnetic  or  any  other  externally  applied  stimuli; 
hence,  interface  can  be  more  practically  described  as  a  region  separating  two  bodies  with 
different  constitutive  response.  Interfaces  have  been  modeled  as  either  continuous  or 
discrete  systems  based  on  the  problem  of  interest  and  the  length  scale.  In  the  case  of 
interfaces  between  two  phases  of  composite  material,  the  primary  interest  is  to  study  the 
load  transfer  between  the  phases  across  the  interface.  For  such  problems,  finite  element 
based  methods  are  adequate  provided  the  general  constitutive  response  of  the  interface 
region  is  known.  Several  such  continuum  methods  have  been  developed;  for  example, 
interface  region  has  been  modeled  as  a  thin  surface  separated  by  springs  [33-35]  and  as 
narrow  continuum  with  graded  mechanical  properties  [36].  Another  method  consists  of 
incorporating  the  interface  within  continuum  description  based  on  cohesive  zone  models 
[37,  38]. 

Interfaces  can  also  be  perceived  as  two-dimensional  defects  in  an  otherwise 
homogenous  atomic  structure.  This  would  be  a  high-energy  region  with  many  processes 
related  to  atomic  structural  changes  and  chemical  reactions  occurring  in  a  narrow  region. 
For  example,  there  is  evidence  of  formation  of  series  of  carbides  in  the  interfaces  of  SiC- 


Ti  metal  matrix  composites  [39],  Atomic  simulation  can  be  employed  to  obtain  the 
constitutive  law  required  for  continuum  formulation  of  interface  incorporating  the 
structural  and  chemical  changes  occurring  at  the  interfaces.  We  now  discuss  the 
representation  of  interfaces  employing  continuum  based  cohesive  zone  models  as  well  as 
discrete  molecular  dynamics  approach. 

Continuum  model  for  interface 

Consider  two  bodies  £2|  and  separated  by  an  inter  facial  region  y  (figure  1 ).  The 
two  bodies  can  be  considered  as  matrix  and  fiber  of  the  composite  material.  In  order  to 
study  the  load  transfer  between  the  bodies  fij  and  Cb.  we  need  to  subject  the  bodies  to  a 
relative  motion  with  respect  to  each  other.  Let  the  system  be  subjected  to  external 
displacement  (or  traction)  as  shown  in  figure  1  (a).  Under  the  action  of  external  load  a 
material  point  initially  at  X  in  either  of  bodies  Qi  or  fi2  moves  to  a  new  location  x.  This 
transformation  can  be  described  by  any  of  the  strain  measures,  for  instance  deformation 
gradient  can  be  expressed  as: 


F  =  — 

dX 


(1) 


The  kinematic  measure  can  be  related  to  its  kinetic  conjugate  through  the  appropriate 
constitutive  relation  as  shown  below 


P  =  CF  (2) 

where,  P  is  the  first  Piola-Kirchoff  stress  and  C  represents  the  linear  or  nonlinear 
stiffness  of  the  material.  Using  this  constitutive  relation  in  the  equilibrium  relation  we 


obtain 


This  equation  can  be  solved  with  appropriate  boundary  conditions  to  completely  describe 
the  bodies  Qi  and  fii  Now,  consider  the  interface  region  y.  A  point  P  on  the  interface  lies 
on  the  surface  described  by  the  unit  normal  N.  If  the  material  deforms  and  fractures  along 
the  interface  as  shown  in  figure  1(b),  two  new  surfaces  are  created  (represented  by  unit 
normals  n  t  and  nj).  The  material  point,  which  was  at  P  in  the  initial  configuration,  now 
corresponds  to  two  points  PI  and  P2.  This  results  in  the  loss  of  one  to  one 
correspondence  between  the  two  (undeformed  and  fractured)  configurations  thus 
violating  the  fundamental  tenets  of  continuity. 

One  of  the  methods  to  counter  this  problem  is  to  use  cohesive  zone  models.  We  assume 
that  the  interface  is  a  different  continuous  material  with  an  independent  constitutive 
relation.  The  constitutive  relation  for  the  interface  material  is  given  as  a  displacement 
dependent  potential  i.e.  the  measure  of  energy  required  to  extend  the  interface  material  by 
a  certain  relative  displacement  8  between  the  two  originally  mating  surfaces.  The 
resulting  tractions  on  the  interface  can  be  expressed  as: 

'/|8|<|SJ>  t  =  o.n  =  /(6)  (4) 

f/|8|>|s,4  t  =  o.n  =  0  (5) 

Now  the  point  P  can  be  considered  as  two  identical  points  PI'  and  P2',  which  after 
deformation  (or  fracture)  move  to  the  new  points  PI  and  P2.  Through  this  method,  we 
can  incorporate  the  interface  deformation  and  fracture  into  the  purview  of  continuum 
mechanics.  This  is  the  advantage  of  CZM  over  ad  hoc  numerical  schemes  such  as 
springs. 

Cohesive  zone  model  was  first  proposed  by  Barrenblatt  [40]  and  subsequently  extended 
to  ductile  and  quasi-brittle  materials  by  Dugdale  [41]  and  Hillerborg  [42]  respectively.  In 


the  recent  years,  CZM  has  been  applied  to  a  number  of  systems  to  solve  numerous 
problems  of  fracture  and  interfaces  in  different  materials.  Different  models  have  been 
used  to  represent  the  traction-displacement  constitutive  relations  and  are  tabulated  by 
Shet  and  Chandra  [38],  The  models  differ  in  the  shape  and  the  points  of  inflexion  i.e.  the 
maximum  load  and  separation  distance.  Chandra  and  coworkers  [37]  emphasize  that  the 
apart  from  maximum  stress  and  separation  distance,  the  shape  of  traction  -displacement 
curve  is  also  an  important  parameter  in  accurately  representing  the  material  behavior. 
Atomistic  model  for  interface 

We  now'  look  at  the  atomic  description  of  the  interface  problem  discussed  earlier  using 
molecular  dynamics.  Molecular  dynamics  involves  numerical  solution  of  New'ton’s  laws 
of  motion  for  a  set  of  atoms.  Here  we  solve  for  the  atomic  degrees  of  freedom  i.e.  the 
positions  and  velocities,  unlike  in  continuum  where  the  solution  is  for  structural  degree  of 
freedom.  Molecular  dynamics  needs  specification  of  the  geometry  of  the  system  (initial 
positions  and  velocities  of  atoms)  and  an  energy  function  that  describes  the  interaction 
between  the  atoms.  For  a  set  of  N  atoms  interacting  through  a  potential  <)>(rjj),  the 
interatomic  force  on  atom  i  due  to  atom  j  (fjj)  can  be  obtained  as: 
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dr 
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The  total  force  experienced  by  atom  i  is  the  sum  of  forces  due  to  the  different  particles  it 
interacts  and  the  forces  external  to  the  system. 


f,=£f«+f0 
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(7) 


The  system  evolves  in  time  and  reaches  equilibrium  and  the  physical  and  structural 
properties  and  the  properties  of  interest  are  calculated  as  statistical  averages. 


Consider  the  same  bodies  £2|  and  fii  separated  by  the  interfacia]  region  y,  but  the 
bodies  Q|.  £2:  and  y  are  now  composed  of  discrete  atoms  interacting  based  on  an 
interatomic  potential.  In  this  case,  the  motion  of  individual  atoms  and  the  forces  acting 
between  them  is  the  basis  for  load  transfer.  Interestingly,  the  concept  of  surface  is  not 
clearly  defined  in  the  discrete  description;  an  envelope  of  all  atoms  bounding  a 
predefined  region  can  be  considered  as  surface.  The  interface  y  can  be  considered  a 
region  bounded  by  two  such  surfaces  (x i  and  X2  in  figure  2). 

In  order  to  displace  the  two  bodies  (fii  and  £2;)  with  respect  to  each  other,  a  lixed 
displacement  (or  force)  is  applied  on  a  set  of  atoms  at  the  end  of  body  £2j  (region  o>  in 
figure  2a)  and  the  effect  is  transmitted  through  the  system.  For  example,  if  the  atom  1 
(see  figure  2b)  is  subject  to  external  displacement  u,  there  is  a  change  in  the  relative 
position  vectors  rjj  as  well  as  the  bond  angles  0|,  with  all  the  j  atoms  it  interacts. 
Correspondingly,  there  the  force  (fij)  experienced  by  the  neighboring  atoms  changes  and 
this  effect  propagates  from  Q]  to  Q2  through  the  interface  region  y. 

The  magnitude  of  the  forces  between  the  atoms  in  the  interface  region  y 
determines  the  extent  to  which  the  external  load  applied  on  body  £2?  is  transferred  to 
body  £2|.  These  forces  are  directly  dependent  on  the  bond  strength  between  the  atoms  in 
the  interface  region.  If  the  bond  strength  of  atoms  in  the  interface  region  is  high,  this 
implies  that  a  small  separation  between  the  atoms  from  equilibrium  bond  length  requires 
a  high  energy  input,  correspondingly  the  forces  between  atoms  that  are  responsible  for 
the  load  transfer  will  be  high.  On  the  contrary,  if  the  bond  strength  in  the  interface  region 
is  low'  (for  example,  Vander wall's  interactions)  the  force  which  is  transferred  to  £2i  form 
£2i  wi II  be  low  resulting  in  a  low  load  transfer. 


In  order  to  obtain  the  traction-displacement  (i-8)  relation  for  CZM  from 
molecular  dynamics,  we  need  to  calculate  the  stress  transmitted  to  the  matrix  (fi|)  for  the 
applied  displacement.  This  is  the  traction  experienced  on  surface  %2  (see  figure  2) 
between  the  interface  and  matrix.  There  are  various  definitions  of  stress  in  MD  and  a 
discussion  regarding  their  applicability  to  CNTs  is  discussed  by  Chandra  and  coworkers 
[43].  Traction  experienced  across  the  surface  can  be  obtained  using  Yip’s  stress  below: 
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Where  a,  P  correspond  to  x,  y  and  z  directions,  i,j  are  atom  indices,  A  is  the  area  of  the 
plane  and  V  is  the  interatomic  potential,  m,  v  and  At  are  the  mass  velocity  and  time 
interval.  The  first  term  corresponds  to  momentum  flux  through  the  area  A  and  the  second 
term  is  total  force  per  unit  area. 

Since  we  are  interested  only  on  the  stress  transmitted  to  the  matrix,  we  need  not 
model  the  matrix  explicitly.  Consider  the  situation  as  shown  in  figure  2c  where  the  matrix 
is  not  explicitly  modeled,  instead  the  atoms  on  the  surface  %2  are  fixed.  The  reaction  force 
experienced  by  the  fixed  atoms  would  be  the  force  transmitted  to  the  matrix  (£1|).  The 
stress  in  equation  (8)  above  can  now  be  calculated  as  the  total  reaction  force  divided  by 
the  corresponding  area. 

The  multiscale  model  proposed  can  be  summarized  as  follows.  Consider  the  situation 
when  Hi  and  fT  can  be  adequately  modeled  by  continuum  methods  connected  by  the 
interfacial  region  y  demanding  atomistic  description  modeled  by  molecular  dynamics. 
The  interfacial  region  can  be  represented  in  continuum  as  a  surface  experiencing  traction 


x  with  displacement  S.  This  traction  can  be  determined  by  molecular  dynamics 
simulations  as  the  force  experienced  on  the  surface  separating  Q2  and  y  per  unit  area.  It 
must  be  noted  that  this  is  a  one  way  coupling  going  from  MD  to  continuum  and  not  vice 
versa.  However,  all  nonlinearities  arising  out  of  geometry,  material  or  boundary 
conditions  at  the  continuum  level  can  be  handled  by  the  model. 

Concerning  the  problem  of  CNT  based  composites  studied  in  this  paper,  the 
interfacial  region  between  CNT  and  matrix  is  modeled  as  hydrocarbon  groups  chemically 
bonded  to  the  CNT.  Chemical  bonding  between  matrix  and  CNTs  is  perceived  to  improve 
strength  immensely;  Frankland  and  coworkers  [44]  have  demonstrated  the  effect  of  cross- 
linking  using  molecular  dynamics  simulation  of  CNTs  in  a  polymer  matrix.  Wagner  and 
coworkers  [11,45]  report  2-2  cycloaddition  reaction  as  a  possible  reason  for  high 
interfacial  strength  in  polyurethane-CNT  composites.  Jia  and  coworkers  [10]  observe  that 
use  of  A1BN  initiator  causes  chemical  bonding  between  PMMA  matrix  and  CNT.  In 
addition,  there  have  been  a  number  of  investigations  of  surface  modification  of  nanotubes 
leading  to  addition  of  various  functional  groups  like  esters,  amides,  hydrocarbons, 
hydrogen,  fluorine  etc  to  CNT  [46,  47],  Based  on  MD  simulations  Namilae  and 
co workers  [48]  have  observed  that  the  mechanical  properties  of  functionalized  CNTs  are 
comparable  to  that  of  defect  free  CNTs.  The  results  of  atomic  simulation  are  now 
discussed. 

3.  Atomic  Simulation  of  CNT  Pullout  Tests 

Evaluation  of  interfacial  properties  in  composite  materials  is  typically  done  using 
single  fiber  pullout  and  pushout  tests.  The  output  of  these  tests  is  the  force  required  for 
debonding  the  fiber  from  matrix  and  the  shear-stress  distribution  along  the  interface.  In 


addition  to  experimental  investigations,  there  have  been  a  number  of  finite  element  based 
numerical  investigations  of  pullout  and  pushout  tests  to  understand  the  thermomechanical 
behavior  of  interfaces.  In  the  case  of  CNT  reinforced  composites,  experimental 

application  of  mechanical  loads  on  single  fibers  is  extremely  difficult  because  of  the 

■ 

nanometer  size  of  fibers.  There  has  been  only  one  experimental  determination  of 
interface  strength  by  pullout  test  using  scanning  probe  microscope  [22 1 .  Conducting  these 
tests  on  a  regular  basis  in  different  materials  for  different  boundary  conditions  is  beyond 
the  current  experimental  capabilities;  hence,  molecular  dynamics  simulations  become  a 
natural  choice  for  studying  the  interface  behavior.  Lord!  and  Yao  [49],  Liao  and 
co workers  [50]  have  simulated  pullout  tests  with  non-bonded  interactions  using 
molecular  statics,  but  there  has  not  been  any  detailed  investigation  in  the  case  of  bonded 
interactions.  We  consider  the  situation  where  matrix  and  nanotube  are  chemically 
bounded  by  hydrocarbons  chains. 

Tersoff-Brenner  potential  [51]  is  used  to  calculate  C-C  and  C-H  chemical 
interactions.  The  potential  has  coordination  dependant  terms,  which  enable  modeling  of 
bond  conjugation  in  carbon,  and  has  been  used  by  numerous  investigators  to  study 
processes  such  as  deformation  of  carbon  nanotubes  [5,  6],  and  various  surface  chemical 
reactions  in  carbon-hydrogen  systems  [52].  Linear  hydrocarbons  modify  the  planar  SP2 
configuration  of  CNT  to  tetrahedral  SP3  bond  with  a  pi  orbital  on  the  CNT  surface.  This 

Q 

change  in  bonding  and  increase  in  bond  length  to  1.54  A  is  accurately  represented  using 
Tersoff-Brenner  potential. 


The  boundary  conditions  applied  for  the  pullout  test  simulation  are  shown 

schematically  in  the  figure  3.  The  corner  atoms  of  the  hydrocarbon  attachments  are  fixed 

0 

indicating  that  they  are  connected  to  matrix  at  those  locations.  Displacement  of  0.05  A  is 

o 

applied  to  the  atoms  at  one  end  of  the  nanotube,  about  15  A  in  length  in  order  to 
simulate  the  effect  of  pullout.  Following  each  displacement,  the  system  was  equilibrated 
for  1500  time  steps  of  0.2  fs  duration.  The  simulations  were  carried  out  until  some  of  the 
hydrocarbon  chains  fail.  Typically,  simulations  lasted  for  500000-800000  time  steps, 
Variations  in  the  length  and  density  (number  per  unit  area)  of  chemical  attachments  were 

investigated  for  a  (10,10)  CNT  122  A  in  length.  The  reaction  force  on  the  fixed  atoms  is 
monitored  throughout  the  simulation  and  is  averaged  over  100  lime  steps  before  applying 
the  next  set  of  displacements. 

Figure  4  shows  typical  force  vs.  displacement  for  any  hydrocarbon  attachment.  The 
force  in  the  figure  is  the  average  reaction  force  experienced  by  the  fixed  atom  in  the 
direction  along  the  length  of  nanotube  and  corresponds  to  shear.  Displacement  is 
calculated  as  the  change  in  the  position  experienced  by  the  atom  that  is  initially  attached 
to  the  CNT.  Though  there  are  statistical  variations  for  different  chemical  attachments,  the 
general  shape  of  the  force  displacement  plot  is  as  shown  in  figure  4. 

The  initial  region  of  curve  is  flat  (parallel  to  displacement  axis)  marked  as  region 
(a).  This  region  corresponds  to  stretching  of  the  hydrocarbon  attachment.  The  flat  region 
shows  that  there  is  minimal  load  transfer  in  this  portion  of  curve  and  is  similar  to  the 
mechanical  analogue  of  a  loose  string  becoming  taut.  The  length  of  ibis  flat  portion  is 
directly  dependant  on  the  length  of  hydrocarbon  chain  i.e.  the  flat  region  is  longer  for 


hydrocarbon  chains  with  four  and  five  carbon  atoms  than  with  two  carbon  atoms.  After 
this  flat  region,  there  is  a  gradual  increase  in  the  reaction  force  corresponding  to  region 
(b)  of  the  curve.  In  this  part,  the  carbon  chain  contributes  to  the  load  transfer.  Though 
there  are  statistical  variations  from  plot  to  plot,  the  typical  force  experienced  in  this 

.  0  , 

portion  of  the  curve  is  about  3eV/A  (4.8  x  10'  pN).  This  value  of  force  is  very  small  but 
it  must  be  noted  that  the  area  on  which  this  force  acts  is  of  the  order  of  angstroms; 
consequently,  the  resulting  shear  stress  is  very  high.  Region  (c)  of  the  plot  as  shown 
consists  of  fluctuations  in  the  reaction  force.  This  is  due  to  an  interesting  behavior  of 
bond  separation  and  rejoining  with  adjacent  atoms;  when  the  separation  occurs  there  is  a 
sudden  drop  in  force  but  this  is  followed  by  rejoining  of  the  hydrocarbon  chain  with 
adjacent  atom  of  the  CNT.  After  the  jagged  region,  there  is  a  sudden  increase  in  reaction 
force  as  shown  in  region  (d)  of  the  figure  followed  by  total  failure.  The  force  at  which  the 
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failure  occurs  is  about  6eV/A  (~10'“  pN).  This  is  the  force  required  to  break  one 
chemical  attachment;  the  overall  force  required  to  break  all  the  chemical  attachments  is 
much  higher  and  is  sum  of  all  individual  reaction  forces. 

The  area  under  the  force-displacement  curve  denotes  the  energy  required  for 
detaching  of  the  hydrocarbon  group  from  CNT.  Calculations  based  on  molecular  statics 
indicate  that  the  energy  required  for  separation  is  of  the  order  of  3  eV.  However,  the  area 
under  the  force  displacement  curve  is  much  greater  and  is  about  20±4  eV  for  various 
attachments.  This  shows  that  the  dynamic  process  of  nanotube  pullout  requires  much 
higher  energy  than  that  predicted  based  on  statics.  Energy  of  3  eV  is  associated  with  each 
time  the  bond  is  broken,  as  this  process  is  repeated  number  of  times  for  each  chemical 
attachment  the  total  energy  consumed  in  the  separation  process  is  much  higher. 


Different  cases  were  considered  with  varying  number  of  chemical  attachments. 
For  the  (10,10)  CNT  studied  here,  there  are  85  chemical  attachments  if  one  chemical 
attachment  is  attached  per  repeat  unit  of  the  CNT  excluding  the  end  of  CNT  that  is 
subject  to  displacement.  Different  number  of  chemical  attachments  85,  45,  20,  10  and  5 
are  attached  at  circumferentially  random  and  cylindrical ly  equidistant  locations  and 
pullout  tests  are  simulated  with  all  the  other  variables  such  as  temperature  and 
displacement  rate  being  constant.  The  general  shape  of  force-displacement  curve  remains 
similar  to  that  shown  in  figure  4.  The  values  of  points  of  inflexion  also  remain  similar 
within  statistical  variation,  but  there  is  a  tendency  for  extended  debonding  and  rebonding 
region  when  the  density  of  chemical  attachments  is  lowered.  This  could  be  ascribed  to  the 
fact  that  the  distribution  of  load  is  uniform  when  fewer  hydrocarbon  chains  are  attached 
to  the  CNT. 

The  debonding-rebonding  behavior  observed  in  pullout  test  is  unique  to  nanoscale 
and  is  generally  not  observed  in  conventional  composites  at  macroscopic  scale. 
Macroscopic  fracture  is  generally  considered  irreversible,  i.e.,  the  energy  required  in 
creation  of  new  surface  during  fracture  is  not  equal  to  energy  required  to  join  two 
surfaces  and  heal  the  damage.  At  atomic  scale,  fracture  is  defined  by  breaking  of 
chemical  bonds.  The  energy  required  to  break  two  bonds  is  exactly  equal  to  the  energy 
released  to  join  the  same  two  atoms  to  form  chemical  bond.  Whether  two  atoms  will  join 
to  form  a  chemical  bond  is  determined  based  on  the  distance  between  the  atoms  and  the 
kinetics  of  the  process.  This  phenomenon  is  not  entirely  new;  for  example,  similar 
observations  have  been  made  during  atomic  simulations  of  stick  slip  friction  and  that  of 
metal  cutting  processes.  Zhang  and  Tanaka  [53]  have  performed  simulations  of  friction 


on  copper  surfaces;  they  observe  variation  of  frictional  force  in  oscillatory  manner  similar 
to  that  of  reaction  force  observed  in  our  simulations.  Further,  the  frictional  force  during 
rolling  of  nanotube  on  graphene  sheet  produces  similar  periodic  variations  [54], 
However,  significance  of  this  phenomenon  is  increased  here  since  it  affects  the  load 
transfer  considerably. 

Traction  displacement  relation  for  nanoscale  interfaces 

We  now  proceed  to  obtain  the  interface  constitutive  relation  based  on  the  above 
discussed  atomic  simulation  results.  An  interesting  aspect  of  using  CZM  to  model  the 
interfaces  is  that,  all  the  micro  mechanical  details  occurring  during  the  load  transfer  and 
failure  are  taken  into  account  in  the  traction  displacement  relation.  For  example,  in 
conventional  fiber  composites  various  micromechanisms  such  as  crack  deflection,  crack 
closure,  fiber  bridging,  inelastic  dissipation  etc  act  simultaneously  and  are  represented  in 
an  average  sense  in  the  extrinsic  and  intrinsic  regions  of  the  t-8  relation.  In  the  case  of 
CNT  based  composites,  various  mechanisms  such  as  chemical  bonding,  bond  breaking 
and  joining,  Vanderwall's  interactions,  nanotube  curvature  etc  may  affect  the  load 
transfer.  By  using  atomic  simulations  to  determine  the  t-5  relation,  we  include  some  of 
these  mechanisms  into  the  continuum  formulations. 

There  is  an  additional  complexity  involved  in  obtaining  cohesive  zone 
information  from  atomic  simulations.  Because  of  discrete  numerical  implementation  of 
finite  element  method,  a  given  discretized  element  experiences  the  constitutive  response 
as  a  single  unit  (albeit  based  on  an  interpolation)  according  to  the  traction  displacement 
response.  The  element  dimensions  can  be  refined  to  the  order  of  few  nanometers,  but  it 
would  still  physically  correspond  to  a  large  number  of  functional  group  attachments.  We 


need  some  way  to  homogenize  the  force  displacement  relations  shown  in  the  preceding 
section  to  obtain  the  traction  displacement  plots  that  are  on  the  scale  of  element.  Simple 
averaging  over  all  the  chemical  chains  will  not  serve  this  purpose  since  the  effect  of 
length  and  density  would  be  ignored.  In  other  words,  if  only  few  chemical  attachments 
near  the  loading  end  take  up  most  of  the  load,  then  the  traction  displacement  plot 
obtained  by  averaging  over  the  entire  length  of  CNT  would  give  an  erroneous 
underestimation. 

Figure  5  shows  variation  of  reaction  forces  for  specific  atoms  in  the  hydrocarbon 
attachments  across  the  length  of  (10,10)  CNT  with  85  attachments.  It  can  be  observed 
that  the  variation  of  force  occurs  in  a  sequential  manner.  We  can  divide  the  typical  force 
displacement  into  various  parts  initial  loading  b,  debonding-rebonding  region  c, 
maximum  loading  region  d,  and  failure  denoted  by  e  (See  figure  4).  Typically,  the  force 
displacement  plots  of  each  chemical  attachment  go  through  these  regions  in  a  self-similar 
manner  until  they  fail.  As  expected  the  chemical  attachments  near  the  loading  region 
reach  the  peak  stress  and  fail  sooner  than  the  chemical  attachments  away  from  the 
loading  end.  Once  the  chemical  attachments  detach  completely  from  the  CNT,  it  fails  to 
take  up  any  more  load;  thus,  the  peak  stress  region  moves  progressively  from  the  loading 
end  to  the  other  end.  Figure  6  schematically  depicts  this  process.  The  various  shades  of 
gray  in  the  figure  6  (b-e)  represent  the  different  stages  of  force-displacement  plot,  from 
no  loading  to  full  loading  and  again  no  load  moving  away  from  the  loading  end.  There  is 
a  gradual  progression  of  the  fractured  region.  For  example.  Figure  5(a)  at  20  and  24  ps 
shows  atoms  near  the  loading  end  taking  up  load.  In  figure  5(b)  the  chemical  attachments 
near  the  loading  end  are  in  debonding-rebonding  phase  while  additional  chemical 


attachments  take  up  load  (region  a).  In  the  figures  5  (c)  and  (d)  the  chemical  attachments 
near  the  loading  end  reach  maximum  load  and  failure  (region  d).  The  maximum  load 
region  moves  away  from  the  loading  end  as  shown  in  figures  5  (e)  and  (f),  At  any  point 
of  time  about  forty  chemical  attachments  only  take  load  and  others  have  either  failed  or 
yet  to  be  loaded.  Based  on  above  observations,  we  have  taken  the  average  of  these  plots 
with  applied  displacement  for  forty  chemical  attachments  and  divided  by  corresponding 
area  to  obtain  the  traction  displacement  plots.  In  the  case  of  lower  number  of  chemical 
attachments,  we  observe  that  all  the  chemical  attachments  take  similar  load  at  any  given 
time  and  the  corresponding  area  is  used  for  the  traction  calculation.  The  traction- 
displacement  plots  for  use  in  cohesive  zone  arc  obtained  by  extrapolating  figures  to 
complete  failure  as  shown  schematically  in  the  figure  7. 

An  interesting  observation  of  this  investigation  is  the  very  high  inlerfacial 
strength  obtained  when  there  is  chemical  bonding  between  matrix  and  CNT.  For  high 
degree  of  chemical  bonding  (85  chemical  attachments)  we  obtain  interfacial  strength  as 
high  as  5  GPa  which  is  about  two  orders  of  magnitude  higher  than  that  typically  observed 

in  conventional  carbon  fiber  reinforced  polymeric  composites  [55].  In  the  case  of  five 
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hydrocarbon  attachments  distributed  over  122  A  we  obtain  an  interfacial  peak  strength 
of  500  MPa.  It  is  interesting  to  note  that  experimental  measurements  of  interfacial 
strength  by  Wagner  et.  al.  [11]  and  by  Cooper  et.  al.  [22]  are  of  the  same  order.  Further, 
Wagner  and  coworkers  indicate  that  this  high  value  might  be  due  to  chemical  interaction 
between  CNTs  and  urethane  matrix.  We  obtain  only  a  value  of  about  50  MPa  of  strength 


when  number  of  chemical  attachments  in  200  A  length  CNT  is  only  three.  These 


observations  indicate  the  strength  of  the  interface  can  be  altered  by  two  orders  of 
magnitude  by  altering  the  density  of  chemical  attachments. 

4.  Cohesive  Zone  Modeling  of  Nanoscale  Interfaces 

We  can  now  apply  the  atomically  (MD)  informed  traction  displacement  plots  in 
finite  element  simulations  of  CNT  reinforced  composites.  In  order  to  demonstrate  the 
effectiveness  of  the  multiscale  model  we  study  the  elastic  behavior  of  the  composite  and 
the  effect  of  interfacial  strength  on  the  stiffness,  it  is  a  common  practice  to  evaluate  the 
mechanical  behavior  of  conventional  composites  using  representative  volume  element 
(RVE).  Because  of  their  ability  to  fill  space,  square  or  hexagonal  RVEs  are  popular  in 
formulations.  Cylindrical  RVEs  using  axisymmetric  formulations  require  less 
computational  effort. 

The  finite  element  mesh  and  the  boundary  conditions  are  as  shown  in  the  figure  8. 
Axisymmetric  four  node  cohesive  zone  elements  with  zero  initial  thickness  in  the 
direction  normal  to  the  interface  are  used  to  model  the  interface  behavior.  The  traction 
displacement  plots  discussed  in  the  previous  section  form  the  constitutive  behavior  of  the 
interface  while  the  matrix  and  the  nanotubes  are  modeled  as  continuum  with  appropriate 
elastic  properties.  Four  node  axisymmetric  elements  are  employed  to  descretize  both 
matrix  and  CNT.  The  models  typically  consist  of  1000  elements  (variation  for  different 
fiber  content)  and  the  interface  consists  of  93  cohesive  zone  elements.  CNT  length  is 
assumed  to  be  450  A  and  diameter  of  20  A.  In  these  computations  CNT  is  considered  as 
solid  cylinder;  this  can  be  justified  by  assuming  a  multiwalled  nanotube  with  innermost 
tube  of  very'  small  radius  (for  e.g.  (5,0)  nanotube).  Both  the  matrix  and  CNT  are  assumed 
continuous  media  exhibiting  homogenous  linear  elastic  isotropic  behavior  with  a  given 


Young’s  modulus  and  Poisson's  ratio.  Young's  modulus  of  1000  GPa  and  Poisson's  ratio 
of  0.3  are  employed  for  the  nanotube.  The  general-purpose  commercial  code  ABAQUS 
[56]  is  employed  to  carry  out  the  analysis.  The  cohesive  element  model  is  input  as  a  user- 
defined  element  subroutine  UEL  into  ABAQUS.  Parametric  studies  for  different 
inlerfacial  strengths  to  understand  the  effect  of  interfacial  strength  on  the  stiffness  and 
load  transfer  of  composite. 

Effect  of  interface  strength 

Two  extreme  cases  of  load  transfer  occur  either  when  the  interfacial  stress  transfer  is 
perfect  or  when  there  is  zero  load  transfer  between  the  matrix  and  fiber.  In  terms  of  CZM, 
these  two  cases  would  correspond  to  infinite  and  zero  peak  traction  respectively.  In  the 
case  of  zero  in  ter  facial  traction,  the  composite  elastic  properties  trivially  become  equal  to 
that  of  matrix  material  as  no  load  transfer  occurs.  The  stiffness  of  the  composite 
progressively  increases  as  the  peak  traction  of  CZM  increases.  For  very  high  value  of 
peak  stress  the  load  transfer  should  be  near  maximum  possible  load  transfer  and 
correspondingly  the  elastic  modulus  of  the  composite  should  be  near  that  of  ideal 
interface.  This  concept  is  illustrated  in  figure  9,  The  figure  shows  variation  of  composite 
elastic  modulus  with  volume  percentage  CNT  for  various  CZMs  described  in  the  earlier 
section.  When  the  matrix  and  CNT  are  connected  by  a  large  number  of  chemical 
attachments  (one  per  repeat  unit  of  CNT),  the  interfacial  strength  is  very  high  and  the 
peak  traction  is  of  the  order  of  5  GPa,  This  high  interfacial  strength  leads  to  stiffness 
values  near  that  of  ideal  composite  with  perfectly  bonded  interface  as  seen  in  figure  9. 

On  the  other  extreme  when  the  peak  traction  is  5  MPa  (this  is  a  fictitious  traction 
displacement  plot  and  does  not  correspond  any  atomic  simulation  result),  the  composite 


stiffness  is  close  to  that  of  matrix  material  (denoted  by  dotted  line  in  figure  9).  The  more 
likely  values  of  inter  facial  peak  traction  are  between  50  and  500  MPa  which  correspond 
to  3  bonds  per  200  A  CNT  and  5  bonds  for  100  A  CNT,  In  these  cases,  stiffness  of 
composite  is  about  80%  and  60%  of  the  maximum  possible  value  (perfect  interface).  This 
emphasizes  the  importance  of  interface  load  transfer  in  determining  the  stiffness  of 
nanotube  reinforced  composites 
Variation  of  fiber  stiffness 

Figure  10  shows  the  variation  of  composite  stiffness  for  various  fibers  (different 
stiffnesses)  and  interface  stiffnesses  in  epoxy  matrix  (E=3.5  GPA).  Different  fiber 
stiffness  values  corresponding  to  those  of  glass  fibers  (E=85  GPa),  carbon  fibers 
(E=300GPa),  SiC  whiskers  (E=400  GPa)  and  CNT  (E=1000GPa)  are  used.  CNT 
reinforced  composites  in  general  exhibit  much  higher  stiffness  than  other  composites 
with  other  reinforcements.  To  provide  a  numerical  perspective  typical  values  of 
interfacial  strength  in  polymeric  composites  are  about  50  MPa,  it  can  be  observed  from 
the  figure  that  CNT  based  composites  with  high  degree  of  chemical  attachments  have  the 
potential  to  reach  stiffness  of  about  300  %  of  that  of  conventional  carbon  fiber  reinforced 
composites.  This  emphasizes  the  potential  properties  possible  by  surface  modification  of 
CNTs. 

Variation  of  matrix  stiffness 

Figure  1 1  shows  the  variation  of  composite  stiffness  with  matrix  stiffness  for  various 
CZMs.  Typical  values  of  elastic  moduli  for  polymeric  systems  range  from  1  GPa  to  10 
GPa,  Four  matrix  elastic  moduli  of  1,  3.5,  5  and  10  GPa  are  considered  in  this  study.  The 
composite  stiffness  increases  with  increase  in  matrix  stiffness  and  decreases  with 


lowering  of  interfacial  strength  as  expected.  The  ratio  of  composite  stiffness  to  matrix 
stiffness  is  higher  for  matrices  with  lower  stiffness.  It  can  be  observed  from  the  figure 
that  low  interfacial  strength  (5  MPa)  results  in  composite  with  stiffness  similar  to  that  of 
matrix  while  composite  with  high  interface  strength  (5GPa)  has  stiffness  close  to  that  of 
perfect  interface.  The  plot  shows  that  high  interface  strength  is  as  critical  as  matrix 
properties  if  not  more. 

5.  Summary 

Carbon  nanotubes,  with  their  unique  mechanical  properties  have  great  potential  as 
reinforcements  in  ultra  strong  and  ultra  stiff  composites.  Weak  interfaces  significantly 
reduce  the  ability  of  CNTs  in  providing  those  high  levels  of  reinforcement.  Chemical 
attachments  to  CNTs  can  provide  high  interfacial  strength.  Summarizing  the  paper: 

■  A  hierarchical  multiscale  model  is  developed  to  study  the  mechanics  of  interfaces 
in  CNT  based  composites.  The  model  employs  atomistically  informed  cohesive 
zone  traction  displacement  relations  to  link  the  atomic  and  macroscopic  scales. 

■  Molecular  dynamics  simulations  of  CNT  pullout  tests  are  performed  to  quantify 
the  interfacial  behavior.  The  interfaces  between  CNT  and  matrix  are  modeled  as 
hydrocarbon  chemical  attachments. 

■  Interfacial  strength  as  high  as  5GPa  can  be  obtained  by  chemically  bonding  the 
matrix  and  CNTs 

■  Interfacial  strength  varies  as  a  function  of  number  of  chemical  attachments  per 
unit  surface  area,  i.e.  density  of  chemical  attachments. 


•  The  multiscale  model  is  employed  to  study  the  effect  of  interface  strength  on  the 
elastic  properties  of  composites.  It  is  found  that  interfaces  significantly  alter  the 
elastic  response  of  CNT  based  composites 
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Figure  1.  Illustration  of  Cohesive  Zone  model  for  interfaces 
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Figure  2*  Atomic  description  of  interface 


Figure  3,  Schematic  of  the  boundary  conditions  applied  in  the  pullout  test  simulation 
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Figure  4.  Reaction-force  vs  displacement  for  a  typical  hydrocarbon  attachment 
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Figure  5.  Variation  of  the  reaction  force  along  the  length  of  (10,10)  CNT  with  eighty  five 
chemical  attachments  at  different  simulation  times 
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Figure  6.  Schematic  showing  the  variation  of  reaction  force  with  time  evolution. 
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Figure  7.  Traction-displacement  plots  obtained  from  atomic  simulation  are  shown  in  (a) 
and  (b).  These  are  extrapolated  to  obtain  the  cohesive  zone  traction  displacement 

relations  as  in  figures  (c)  and  (d) 


Interface  elements 


Figure  8.  Typical  mesh  used  in  the  finite  element  simulations 
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Figure  9.  Variation  of  composite  Young’s  modulus  with  volume  %  CNT  for  different 

interfacial  strengths. 
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Figure  10.  Variation  of  composite  Young’s  modulus  with  fiber  stiffness  for  different 

interfacial  strengths. 
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Figure]  I .  Variation  of  composite  Young’s  modulus  with  matrix  stiffness  for  different 

interfacial  strengths. 


